Multiscale Entanglement Renormalization Ansatz for Kondo Problem 
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We derive the multiscale entanglement renormalization ansatz (MERA) for the single impuity 
Kondo model. We find two types of hidden quantum entanglement: one comes from a finite- 
temperature effect on the geometry of the MERA network, and the other represents screening of the 
impurity by conduction electrons. As the latter starts to dominate the electronic state, the Kondo 
physics emerges. The present result is a simple and beautiful example of a holographic dual of a 
boundary conformal field theory. 
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The multiscale entanglement renormalization ansatz 
(MERA) offers the best variational wave function for 
quantum critical systems [l|. The MERA is also recog- 
nized to be a discrete version of the anti-de Sitter space / 
conformal field theory (AdS/CFT) correspondence that 
has been a hot topic in superstring theory @, [![. Thus, 
deeper understanding of its functionalities similar to 
those of AdS/CFT and constructing efficient analytical 
or numerical methods for real applications to condensed 
matter physics are necessary. 

In this paper, we apply MERA to the single impu- 
rity Kondo model [U, y5[ . The MERA network repre- 
sents a holographic space where quantum entanglement 
propagates. By introducing a finite-temperature MERA 
network, we can clearly see the entanglement struc- 
ture between the impurity and conduction electrons @. 
Then, the Kondo temperature Tk can be obtained from 
the MERA network itself without assuming the specific 
Kondo interaction. Therefore, Tk is a universal pa- 
rameter determined from the global geometric structure 
of the holographic space. Furthermore, in connection 
with AdS/CFT, we argue that the Kondo physics is one 
of condensed-matter realization of black hole thermody- 
namics. Since the Kondo Hamiltonian is the most typical 
correlated electron system and its critical behaviors have 
been well known, some string theorists have recently at- 
tacked this problem by means of AdS/CFT How- 
ever, the theoretical treatments are in some sense far 
from condensed-matter terminologies. Our work trans- 
lates them into much simpler MERA language. 

We start with the impurity Kondo model in spatially 
one dimension which is a good example to study bound- 
ary critical phenomena and the local Fermi liquid theory. 
The Hamiltonian is given by 



with spin S, Uj is electron hopping, and J is Kondo cou- 
pling. 




FIG. 1: The MERA network in the impurity Kondo model. 
Filled dots, diamonds, and filled triangles are the original sites 
in a quantum ID critical system, disentangler tensors, and 
isometries (except for the top tensor), respectivery. An open 
circle represents the impurity site, and open triangles (shaded 
area) are boundary tensors. Note that the vertical direction 
denotes renormalization flow. 

We consider a variational ansatz appropriate for the 
Kondo Hamiltonian. Because of the criticality, the ansatz 
has a MERA network with some modification near the 
impurity site. Figure [T] is the graphical representation 
of the MERA network. The impurity is located at the 
edge of half-infinite chain, but the impurity can inter- 
act with the conduction electrons at any site through the 
holographic branches. Thus, we believe that this config- 
uration is a good starting point. The variational ansatz 
in the ground state is represented as 
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where Cj (cj ) is spinor representation of an electron anni- 
hilation (creation) operator at site i, S is a local impurity 
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where s is degree of freedom of the impurity, and {rrij} 
are degrees of freedom of the conduction electrons. The 
coefficient (j s,mim2 "' mL can be decomposed into a prod- 
uct of three tensors 



C 



s,mim2 ■ ■■rriL 



E 



a OL\OL>2, ■■■CKjv OL\OL2 ■ ■ "CKjV 



2 



where the number of layers is given by N — log 2 L — 1, 



T: 



is the top tensor, $™.Vj^ Q 



is the MERA network 
1, 2, TV) and a, and the 
boundary tensor Ba^a 2 ---a N (like the 'skin' of an onion 
<&) is defined by the matrix product state (MPS) form 
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where W is the isometry tensor. The indices {ctj} rep- 
resent how the impurity interacts with the conduction 
electrons. Here we take ctj = 1,2, ...,Mj. Note that the 
boundary tensor connects with the MERA network of 
the conduction electrons only through the disentangler 
tensors and the top tensor as shown in Fig. [T] In addi- 
tion, the boundary tensor does not contain disentanglers. 
In the continuous limit, the distance between the near- 
est neighbor isometries at the boundary is different from 
that in the bulk, since all of the isometries and the disen- 
tanglers are the same lattice points. This will lead to the 
difference between the bulk and the boundary scaling di- 
mensions. Thus, our MERA network naturally describes 
the boundary CFT feature [HjJ]. 

Let us move to finite-temperature cases. In our recent 
paper associated with finite-temeperature MERA, it is 
possible to give a gravitational interpretation to the ge- 
ometry of the holographic space [6| . As shown below, the 
MERA network is doubled, and two copies are combined 
with each other in IR region of the network. Then, we 
have found that the interface of the combined networks 
can be regarded as an event horizon. Let us imagine 
that this statement is also meaningful for the interface 

(we suppose that 
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represents the impurity degree of freedom). 



This means that the impurity behaves as a black hole. 
Physically, this would be reasonable, since the screen- 
ing of the impurity by the conduction electrons occurs, 
and some information associated with the spin degrees of 
freedom is lost due to their singlet formation. Since the 
singlet is the maximally entangled state, the maximally 
entanglement entropy is actually equivalent to the black 
hole entropy pl-[ia|. 

At finite temperatures, the wave function is described 
by using thermofield dynamics [l9l - l24j . As already exam- 
ined in the previous paper Q, the thermal state is given 
by the MERA - tilde MERA combined network: 
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and C^'™ im2 --- mL is truncated MERA network [2J, [2€ 
defined by 
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where N is the number of layers along the vertical direc- 
tion. It is noted that the index a is put onto the MERA 
network for the conduction electrons. 




FIG. 2: Dual MERA network. A red wavy line represents an 
event horizon arizing from finite-temperature effect, while a 
dashed line represents a horizon associated with entanglement 
between the impurity and the conduction electrons. 




FIG. 3: Spin correlation in the holographic space. When the 
horizon coming from the impurity is terminated by the wavy 
horizon, the long range spin correlation is lost. 



Figure[2]is the graphical representation of the network. 
We find that the indices a and u> represent degrees of 
freedom of the event horizon arizing from finite temper- 
ature effect (we may neglect the effect of uj on the en- 
tanglement entropy for large L). On the other hand, 
the a parameters represent degrees of freedom for entan- 
glement between the local impurity and the conduction 
electrons. Therefore, two types of different horizons ap- 
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pear. Strictly saying, the latter is artificial one, since 
the black hole should contribute to thermalization of the 
system that only comes from the former one as shown 
later. 

When the latter dominates the electronic states, the 
impurity spin can be efficiently screened by the conduc- 
tion electrons. This is holographic representation of the 
Kondo effect. Figure [3] shows the spin correlation in the 
holographic space. At zero temperature (Fig. [IJ , the cor- 
relation is long-ranged. As temperature increases, some 
upper layers of the MERA network start to truncate. 
Then, the wavy horizon terminates the impurity horizon, 
and the long-range spin correlation is lost. The interest- 
ing thing is that we do not explicitely treat the Hamil- 
tonian. Of course, our MERA network has information 
about the location of the impurity and the singlet forma- 
tion among the impurity and the conduction electrons. 
However, the information is indirect. We think that even 
such loose conditions are enough for describing critical 
behavior of this system. 

In the following, we determine the Kondo temperature 
Tk from the MERA network. For this purpose, we calcu- 
late the black hole entropy by based on the holographic 
theory [f| [53, HI] . The entropy is given by the logarithm 
of the total degrees of freedom at the two horizons \ as 
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Here, our original discrete network in Fig. [T]has the fol- 
lowing metric in the continuous limit (rj = 2 in a case of 
binary MERA network): 



ds 2 = (drhirj) + {^ T dx) . 
When we dehne the variable z as 

z = if , 

we obtain the standard notation of the AdS space 
I 



ds 2 = — [dz 2 + dx 2 ^ 



In the present case, x is represented as 
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where N(= tq) is the bond number at the interface be- 
tween B and $ (or the layer number of the MERA net- 
work, since the bond number is 1 per layer), and A is 
the bond number at the MERA - tilde MERA interface. 
When the impurity is disentangled from the conduction 
electrons (Mj — 1 for all j), the previous result is ob- 
tained [6,]. Here, m is the rank of the tensor at the latter 
interface. The area A is given by 
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and for simplitity we take 77 = e, although our figures are 
binary MERA networks. The paramter m is the tensor 
dimension at the horizon. Then, we find 



N 



Sbh = — lnm + y^lnMj. 

6 3=1 



(13) 



By changing the variable To to zq according to Eq. (0), 
Sbh is finally represented as 



S BH = -^flnm + Alnm + ^ln (—1 J (14) 
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where the bulk part Sbulk arid the boundary part 
Sboundary are respectively defined by 



In m, 



(15) 



Sbulk = l*(fe LM) ) 
Sboundary = In 2 In m + In ( — J ■ (16) 

(We may have some umbiguity for the determination of 
the boundary part due to the discretization of our space- 
time.) 

Equation (fTS"]) should be consistent with the bulk en- 
tanglement entropy evaluated by the CFT living in the 
UV boundary of the MERA network [13, Hi]. The en- 
tropy is actually given by 



S 



EE 



— In ( — sinh ( — 
6 



7T 



c (fi firL 



(17) 



for small T with a finite temperature resolution L _1 [29], 
l30j . Here, c is the central charge, and the UV cutoff is 
taken to be unity. Identifying Eq. (IT51) with Eq. (|T7|) . we 
find the temperature of this system as 



k B T = 
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Equation 



comes from the g function of the bound- 
ary CFT |31| . It is clear that we have the following 
relations at the horizon when the three indices of the 
boundary isometry have the same dimensions: 



M 



N 



(19) 



Thus, the result indicates that the boundary entropy 
ln(Mf/m) at layer j decreases under the boundary RG 
flow. This seems to be consistent with the conjectured g 
theorem. Since the IR behavior of the Kondo model is 
described by the local Fermi liquid [32[ , the effect of the 
impurity on the network is smeared out in the IR region. 
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The screening of the impurity spin becomes dominant, 
when the first term in Eq. (fLU)) is greater than the sec- 
ond term. The boundary is determined by solving the 
following equation for Nk 

L Nk 

In m = 5>M„ (20) 

j=i 

and the Kondo temperature is given by Eq. (fl8]) with use 
oiN K {=T K ): 

k B T K = — = ~e- TK . (21) 

TTZk 7T 

Let us compare this result with the well-known Kondo 
temperature UbTk ~ exp (— 1/Jp) where p is the density 
of states. It is hard to solve Eq. (I2TH) . but it's easy to 
find that tk increases with L. Thus, tk increases with 
decreasing the impurity density and p. Actually, this 
feature is not contradict to the exact Tk- 




(b) 



FIG. 4: Holographic representation of the RKKY interac- 
tion: (a) low termerature case, (b) high temperature case. In 
Fig.(b), two impurities can not interact with each other. Red 
lines are bending in these figures, but they are geodesic lines. 

The consistency between the main result Eq. (|2"T|) and 
the exact Tk shows the powerfulness of the MERA anal- 
ysis. The present method can be easily applied to other 



systems. For instance, let us consider the competition be- 
tween the Kondo singlet formation and the Ruderman- 
Kittel-Kasuya-Yoshida (RKKY) internaction in a two- 
impurities case 33]. We assume that the impurities are 
located at both edges in the ID chain. Then, the RKKY 
interaction occurs through the holographic tree. When 
we consider finite temperature cases, some upper layers 
of the MERA network are terminated. Then, the dis- 
tance between two impurities that can interact with each 
other is restricted. As shown in Fig.QJb), two impurities 
can not contanct with each other at high temperature 
cases, since the two geodesic curves though the conduc- 
tion electron is too short to combine these impurities. 
At the same time, the partial screening of each impurity 
is originated only from near-neighbor electrons. Even in 
such a case, the Kondo coupling occurs for large J, since 
J is originally local interaction. On the other hand, at 
low temperatues and for relatively small J values, the 
RKKY interaction wins. As shown in Fig. Ufa), the free 
energy for RKKY is lower than that for Kondo interac- 
tion, since the thermal entropy (or geodesic distance) is 
larger in the former case. Our view seems to be consistent 
with the scaling relations Tk ~ e _1 / J and Trkky ~ J 2 
for weak J. 

Summarizing, we have found two types of entangle- 
ment in the MERA network of the impurity Kondo 
model: one comes from finite-temperature effects on the 
network geometry, and the other represents the screen- 
ing of the impurity by the conduction electrons. As the 
latter starts to dominate the electronic state, the Kondo 
physics emerges. We determine the Kondo temperature 
Tk in a gravitational viewpoint where Eq. f| 12[) plays an 
important role on connecting the network with the di- 
crete AdS space. Since Tk is essentially an universal 
parameter, this temperature scale automatically appears 
from the network geometry itself. Quite recently, it has 
been discussed that the one-particle Schrodinger equa- 
tion for the Wilson Hamiltonian is similar to a bound- 
ary field theory in the background AdS3 spacetime [3~ij |. 
These recent results illustrate close relationship between 
the renormalization group approach in condensed matter 
physics and the AdS/CFT correspondence. 
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